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ABSTRACT 

We carry out plasma diagnostic analyses for 123 planetary nebulae (PNe) and 42 
H II regions using the N II and O II optical recombination lines (ORLs). New effective 
recombination coefficients for the N n and O II optical recombination spectra are used. 
These data were calculated under the intermediate coupling scheme for a number of 
electron temperature (T e ) and density (N c ) cases. We used a new method to determine 
the T e 's and N e 's for the nebular sample, combining the ORLs with the most reliable 
measurements for each ion and the predicted intensities that are based on the new 
atomic data. Uncertainties of the derived T c and N c are estimated for each object. 

The diagnostic results from heavy element ORLs show reasonable agreement with 
previous calculations in the literature. We compare the electron temperatures derived 
from the N n and O n ORLs, T (ORLs), and those from the collisionally excited 
lines (CELs), T (CELs), as well as the hydrogen Balmer jump, T C (H I BJ), especially 
for the PNe with large abundance discrepancies. Temperatures from He I recombi- 
nation lines, T c (He i), are also used for comparison if available. For all the objects 
included in our sample, T c (ORLs) are lower than T C (H I BJ), which are in turn sys- 
tematically lower than T c (CELs). Nebulae with T c (He i) available show the relation 
T c (ORLs) < T c (He i) T C (H I BJ) < T c (CELs), wh ich is cons i stent with predictions 
from the bi-abundance nebular model postulated bv lLiu et al.l (|2000f ). 

Key words: atomic data - planetary nebulae: general - H II regions 



1 INTRODUCTION 

Photoionized nebulae, such as planetary nebulae (PNe) and 
H II regions, provide much of our knowledge of elemental 
abundances in the Milky Way and other galaxies. Accurate 
measurements of electron temperatures (T e ) and densities 
(N e ) are essential for reliable determination of elemental 
abundances. Until recently, the principal means of deter- 
mining elemental abundances in nebulae has been from the 
measurement of collisionally excited lines (CELs). The emis- 
sivities of these lines relative to a hydrogen line are very 
sensitive to T c under typical physical conditions of pho- 
toionized nebulae, and thus are much affected by the errors 
in electron temperatures. These electron temperatures are 
in turn deri ved from the traditional me thod based on the 
CEL ratios l|Osterbrock fc Ferland 2006), such as the [O III] 
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A4363/(A4959 + A5007) and [N n] A5755/(A6548 + A6584) 
nebular-to-auroral line ratios. 



An alternative method of nebular abundance determi- 
nations is to ratio the intensity of an optical recombination 
line (ORL) of helium or a heavy element with that of hy- 
drogen. Unlike CELs, such as the [O ill] and [N n] nebular 
lines, whose emissivities relative to a hydrogen recombina- 
tion line increase exponentially with T a , the emissivities of 
heavy element ORLs relative to a hydrogen recombination 
line change weakly with both T c and N c , as the emissivi- 
ties of recombination lines have only a similar, power-law 
dependence on T B and N e , apart from the parentage effects 
to be discussed in Section 2.2, and are essentially indepen- 
dent of N e under typical nebular physical conditions (e.g. 
iLiul l2006 a, b). Consequently, the method based on ORLs 
is much less affected by temperature measurement errors, 
and the results in principal should be more conclusive. It is 
possible to use the intensity ratio of two ORLs originating 
fro m the recombination of dif ferent ion parents to estimate 
N e l|Fang. Storey fe Liull201lh . since the density sensitivity 
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of the parent populations is partially reflected in the resul- 
tant ORL relative intensities. 

While emissivities of heavy element ORLs have only a 
relatively weak, power-law dependence on T e , this depen- 
dence varies for lines originating from levels of different or- 
bital angular momentum quantum number I. Therefore, the 
relative intensities of ORLs can be used to derive electron 
temperatur e, provided that very accurate measurements ca n 
be secured (|Liull2003l ; [Liu" et al. 1 12004 iTsamis et all 12004 ). 
Since the sensitivity of the ORL ratios to electron temper- 
ature and density is very weak, in order to obtain T e 's and 
N c 's of the nebular regions where ORLs arise, one needs to 
acquire high-precision measurements of these ORLs. How- 
ever, given the low nebular abundances of heavy elements 
(~10" 4 -f0 -3 or even lower relative to hydrogen) and the rel- 
atively long time scales for a heavy element ion to recombine 
with an electron under the physical conditions of gaseous 
nebulae ijOsterbrock fc Ferlandl 120061 ). the heavy element 
ORLs are very weak compared to the CELs and low-order 
hydrogen lines. Therefore, high-precision measurements of 
those ORLs are required. 

Although several deep spectroscopic surveys have 
been carried out during the past decade for several 
dozen Galactic disk an d Bulge PNe |Tsamis et al. [2003b, 



20041; iLiu et al.l |2004|; IRobertson-Tessi fc Garnettl 12005 



Wesson et al. 1 120051 ; IWang fc Liu] I2007T I and for a number 



of G al actic a nd extragalac t ic H h regions jEsteban et al.l 



20021. |2004 ITsamis et all l2003al; iPeimbert et all |2004 



Garcfa-Roias et al.f2 004. 2005, 200<t), most of the ORL mea- 



surements are still not accurate enough for nebular anal- 
ysis, due to either relatively low signal-to-noise ratios or 
line blending. Currently only a limited number of objects 
(mainly nearby Galactic PNe) h ave deep enough spec t ra for 
recombination line analys i s fe.g.lLiu et al]|l995l. [20001. |200H; 
Esteban et"al"1ll999|. l2004 ISharpee et"afll2004 IZhang et all 
2005al ; lFang fc Liulhoill ). 



In nebular astrophysics, there has been a long-standing 
dichotomy in abundance determinations and plasma diag- 
nostics. For a given PN, the abundances of heavy-element 
ions (C 2+ , N 2+ , 2+ , and Ne 2+ ) relative to hydrogen derived 
from the nebular ORLs are all higher than the correspond- 
ing abundance values derived from their CEL counterparts. 
The electron temperatures derived from the H I Balmer 
jump are also systematically lower than those derived from 
the CELs. A number of mec hanisms have been proposed 
to explain this dichot omy (e.g. |Peimberd[l9"o7l ; iRubinll 19891 ; 
IViegas fc cieg3ll994h . but have failed to provide a consis- 
tent interpretation of all observations, especially for those 
PNe with dramatically large abundance discrepancies (e.g. 
>10, 20) . A b i-abundance nebular model postulated by 
ILiu et ail (2000) provides a more natural explanation of this 
dichotomy. In this model, the ORLs heavy element ions 
arise mainly from "cold" H-deficient inclusions, while the 
strong CELs are emitted predominantly from the warmer 
ambient ionised gas of "normal" (~ solar val ues) chemical 



compo s ition. Deep spectroscop i c surveys (e.g . Tsami s et al 
2003b, 120041; ILiu et al.l 12004; IRobertson-Tessi fc Garnett 



20051 ; I Wesson et al.|[2005l ; I Wang fc Liull200it ) and ORL anal- 
ysis of individual nebulae in the past decade (e.g. ILiu et all 
Il995l , |2000| . 120011 1200681 ) have yielded strong evidence for 
the existence of such a " cold" component (Rece nt reviews 
on this topic are given bv lLiull2003l , l2006bl . I20TII '). More re- 



cently, iNicholls et all (|2012f ) proposed that a K-distribution 
for electron energies, which is a departure from the conven- 
tional assumption of a Maxwell-Boltzmann equilibrium en- 
ergy distribution, can explain the dichotomy in both H 11 re- 
gions and PNe. 

Despite the paucity of high-quality observational data 
of the heavy element ORLs as well as the lack of the 
atomic data that are adequate enough for nebular analy- 
sis, plasma diagnostics have been carried out for a number 
of PNe and preliminary results have been obtained. Using 
the ra d iative recomb i nation coefficients of IPequignot et al] 
l|l99ll ). Frsamis et al] (|2004T ) derived electron temperatures 
using the O 11 ORLs for 6 Galactic PNe, and found an 
average T c (0 11 ORLs) value of 2920 ± 2690 K. ILiu et al] 



(2004) derived T c (0 n ORLs) fo r 18 PNe, and foun d an 
average value of 4910 ± 4060 K. IWesson et all (120051 ) de- 
rived electron temperatures for a dozen Norther Galactic 
PNe using the O n ORL ratios A4075/A4089, A4649/A4089 
A4072 /A4089, A4414/A4089 and V1/A4089. IWang fc Liul 
|2007l ) determined T c (0 n ORLs) for 11 Galactic PNe using 
the O n A4089/A4649 line ratio, and gave an average value 
of 4370 ± 5760 K. However, many of the recombination line 
intensities used for diagnostics were not corrected for line 
blending, and the derived electron temperatures could be 
problematic. 

In this paper, we present a new approach to plasma 
diagnostics using the N 11 and O 11 ORLs for a sample of 
PNe and H 11 regions, which are selected from previous op- 
tical recombination line surveys and individual-object stud- 
ies. The most recent effective recombination coefficients are 
utilised for the analysis. For each ion, we use one set of 
detected ORLs with the most reliable measurements to con- 
strain T and iV c simultaneously, by comparing the observed 
intensities with the predicted ones. The primary purpose 
of the current paper is to show that the plasma diagnostic 
method employed here is applicable to all emission nebulae, 
provided that more and more accurate measurements of the 
N 11 and O 11 ORLs are available in the future. 



2 OPTICAL RECOMBINATION LINE 
DIAGNOSTICS 

In this section, we compare the intensities of the observed 
N 11 and O 11 lines with theoretical predictions over a wide 
range of electron temperatures and densities. 

2.1 Planetary nebula and H 11 region sample 

During the past two decades, several deep optical spectro- 
scopic surveys have been carried out and published for sev- 
eral dozens of Galactic disk, Bulge and Halo PNe and for a 
number of Galactic and extragalactic H 11 regions. These al- 
low for detailed nebular plasma diagnostics and abundance 
analyses using the relatively weak hydrogen and helium re- 
combination lines/continua and ORLs from heavy element 
ions. In total, over 100 PNe and 40 H 11 regions have been 
studied using ORLs. The current paper makes use of the 
samples from the literature listed in Tables 1 and 2. 

High spectral and spatial resolution spectroscopy of 
PNe has been published in the literature for detailed analy- 
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Figure 1. Theoretical fractional intensities of O II Multiplet Ml: 2s 2 2p 2 3p 4 D° - 2s 2 2p 2 3s 4 P as a function of electron density. The 
numbers in the brackets (J2 — Ji) following the line labels are the total angular momentum quantum numbers from the upper to the 
lower levels. Transitions from the upper levels with the same angular momentum quantum number J2 are represented by the curves 
with the same colour and line type. Four temperature cases, log T e [K] = 2.5, 3.0, 3.5 and 4.0, are presented. The figure is based on 
the unpublished calculations of Storey. 




Figure 2. Loci of the O II recombination line ratios I(A4649)/I(A4662) and I(A4649)/I(A4089) for different T e 's and N c 's. The figure 
is based on the unpublished calculations of Storey. 
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Figure 3. The O II A4649/A4662 recombination line ratio as a 
function of electron density. Different curves represent different 
temperature cases. The plot is based on the unpublished calcula- 
tions of P. J. Storey. 

of the physical condition an d elemental abunda nce dis- 
tributions across the nebulae. iTsamis et al.l (|2008f ) carried 
out a dedicated study of three Galactic PNe (NGC5882, 
NGC6153 and NGC7009) by means of optical inte- 
gral field spectroscopy using the Fibre Large Array 
Multi Element Spectrograph on the Very Large Telescope 
(VLT/FLAMES). The ratio of the abundances derived from 
ORLs and CELs, which is called the abundance discrepancy 
factor (ADF), was studied across the nebulae for doubly- 
ionised oxygen, 2+ . Very small val ues of the temperature 
fluctuation parameter t 2 (defined by |Peimbertlll967l . Il97lf ) 
in the pl ane of the sky were found in all three objects. Most 
recently, iMesa-Delgado et al.l l|2012f ) published results from 
integral field spectroscopy of a field located near the Trapez- 
ium Cluster in the Orion Nebula using the Potsdam Multi- 
Aperture Spectrophotometer (PMAS). Detailed studies of 
one of the three most prominent protoplanetary disks (pro- 
plyds) show that the ADF of 2+ is close to one in the pro- 
plyd, while the background emission still yields the typical 
ADF(0 2+ ) observed in the Orion nebula. 



2.2 Atomic data 

Since the 1990s, it has been possible to obtain reliable mea- 
surements of the faint ORLs emitted by heavy element ions 
for bright nebulae. However, atomic data are necessary to 
analyse these spectral features, in particular the effective 
recombination coefficients. ORL ratios from states of dif- 
ferent orbital angular momenta do show some temperature 
dependence and therefore can be used to measure the a ver- 
age t emperature under which the lines are emitted (e.g. iLiul 
l2003h . 

Relative populations of the fine-structure levels of the 
ground term of a recombining ion, such as N 2+ in the case 
of N n and 2+ in the case of O n, deviate from the Boltz- 
mann distribution and vary as a function of electron density 



[ii- 
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)1; 



in low-density nebulae. This variation in the level population 
is reflected in the emissivities, and thus in the effective re- 
combination coefficients of the recombination lines that are 
formed from different parent levels. In essence, the intensity 
ratio of two such lines can be used t o determine electron 
density fe.g. lFang, Storey fe Liull201ll V 

High-precision quantum mechanical calcula- 
tions of the effective recom bination coefficients for 
N II (|Fang. Storey fc Liu! 1201 if ), and O II lines (Storey, 
unpublished) have been completed. Both calculations im- 
prove t hose from previous w ork (e . g. iNussbaumer fc Store 
1 1984 lEscalante fc Victoif Il990l; iPecmignot et all Il991l 
I Storevi 11994 iLiu et al.l 1 19951 ; IKisielius etluT Il998f 
iKisielius fc Storevf l2002). The strongest and best-observed 
lines for N II so far have been those from multiplets V3 
3p 3 D - 3s 3 P° and V39a,b 4f G[7/2,9/2] - 3d 3 F°. For O II, 
they are from multiplets VI 3p 4 D° - 3s 4 P, V10 3d 4 F - 
3p 4 D°, and V48a,b 4f G[5,4]° - 3d 4 F for O II. 

Fig. 1 shows the theoretical fractional intensities of 
the fine-structure components of the O II Multiplet VI 
2p 2 3p 4 D° - 2p 2 3s 4 P as a function of electron den- 
sity at four temperature cases. The same can be seen 
for t he N H multiplet V3 2p3 p 3 D - 2p3s 3 P° in Fig. 5 
from iFang. Storey fc Liuf (120111 ) . Fig. 2 shows the loci of 
the O II recombination line ratios I(A4649)/I(A4662) and 
I(A4649)/I(A4089) for different T c 's and jV e 's. The same can 
be seen for N II recombination line ratios I (A5679)/I(A5666) 
and I (A5679)/I(A4041) in Fig. 8 from IFang. Storey fc Liu] 
(2011). Using the accurate measurements of the N II and 
O II line ratios, we can determine T c 's and JV e 's simultane- 
ously from the loci. Since different ORL ratios of N II or 
O II have different sensitivity on T e and N e , this will affect 
the reliability of the results. 

Given the different temperature-dependence of diag- 
nostic lines, individual diagnostic line ratios are expected 
to yield differing results if the nebulae are not isothermal. 
While the small values of t 2 appear to suggest there are 
no large temperature variations where the CELs are emit- 
ted, the discrepancies between the CEL and ORL diagnos- 
tic results definitely exist, thereby suggesting that the gas 
is not isothermal. While the authors believe that the bi- 
abundance nebular model is a satisfactory explanation of the 
abundance (and temperature) discrepancies, an alternative 
explanation involvi ng non-equilibrium e l ectron energies has 
also b een proposed (Nicholls ct al. 20121; lOwocki fc Scudded 
Il983h . According to lNicholls et al.1 l|2012f ). physical processes 
such as magnetic reconnection, injection of high-energy elec- 
trons through photoionization by a "hard" photon spectrum, 
may generate such non-equilibrium electrons. The role of 
Shockwaves in generating the discrepancies also needs to be 
explored. The choice of abundance discrepancy model, how- 
ever, does not impact on the validity of the ORL diagnostics 
proposed here. 

The traditional method, which uses one line ratio to 
determine an electron temperature or a density, has the 
disadvantage that the density or temperature is undefined. 
Although a temperature-sensitive recombination line ra- 
tio is usually quite insensitive to electron density, some 
density-sensitive line ratios, to a noticeable extent, may 
vary with electron temperature. As an example, Fig. 3 shows 
the O II A4649/A4662 ratio as a function of electron den- 
sity. From the Figure one can see that at low densities 
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(log T c < 3.6), the electron density yielded by the observed 
line ratio may differ by as much as 0.6 dex when the logarith- 
mic temperature log T c increases from 2.6 to 3.6. Instead of 
using a single line ratio, we present our new method using 
one set of N II or O II ORLs to constrain the electron tem- 
perature and density simultaneously. Our method is based 
on the fact that the intensity ratio of two of these selected 
ORLs is temperature-sensitive, while the ratio of the other 
(two) lines of this group should be density-sensitive. 



2.3 Theoretical intensities 

For the N II effecti ve recombination coefficients 
|Fang. Storey fc Liul l201ll ). log T c [K] covers a range 
of 2.1 to 4.3, with an incremental step of 0.1, and 
log N e [cm -3 ] spans from 2.0 to 6.0, with an incremental 
step of 0.1. For the O II effective recombination coefficients 
(Storey, unpublished), log T c [K] covers a range of 2.6 to 4.2, 
with an incremental step of 0.2, and log A c [cm -3 ] spans 
from 2.0 to 5.0, with an incremental step of 0.2. Both 
theoretical and observed intensities are normalised such 
that H/3 = 100. 

For most objects, the other O II multiplets are not as 
strong as those considered in Figs. 1 and Fig. 2, and are not 
suitable to constrain T c and N e , due to data quality. Figs. 4 
and 5 show the emissivities of N II V3 and V39 ORLs (V3 
A5666, V3 A5679, V39b A4041 and V39a A4035) and O n VI 
and V48 ORLs (VI A4649, VI A4662, V48a A4089 and V48c 
A4087), respectively, as a function of T e and A e , where the 
emissivity is derived from the effective recombination coef- 
ficients as follows: 



£ (A) = a cff (A) x (1) 

with h as Planck's constant and c as the speed of light, all 
in cgs units. Due to coarse grid resolution, we performed a 
bilinear interpolation on the effective recombination coeffi- 
cients of the O II ORLs provided by Storey (unpublished). 
Bilinear interpolation simply extends linear interpolation for 
functions of two variables, i.e. a e ff (T e , A c ). 

Based on the calculations of the effective recombination 
coefficients for the N II and O II lines, theoretical intensity 
for an N II or O II ORL relative to H/3 can be calculated as 
a function of T c and A c as follows: 



Ipred (A) 



I (A) _ Q cff (A) 4861 X+ 
I (H/3) ~ a e3 (H/3) A H+ 



x 100, 



(2) 



where I pre d (T C ,N C ) is the theoretical predicted intensity 
of the transition A, normalised such that /(H/3) = 100. 
a e ff(A; T e ,N e ) is the effective recombination coefficient for 
the transition A, and a e g(H/3; T e ,N e ) is the effective re- 
combination coefficient for H/3, which is adop t ed fro m the 
hydrogenic calculations of IStorev fc Hummer] fl995). The 
ionic abundances relative to hydrogen derived from ORLs, 
X + /H + , can be found in the literature for each nebula. 

If the spectral resolution for the nebula is relatively 
poor, for instance FWHM < 1.5 A, then some blending may 
occur. The theoretical intensities from the blended lines are 
then added together in the following manner: 



i=l 



Lpred 



(A*) = 



a off (A t ) N, 4861 X+ 
Q cff (H/3) H+ 



x 100, 



(3) 



where i goes from 1 to the number of the blended lines. 

2.4 Plasma diagnostics 

2.4-1 Least-squares fit 

For each PN and H II region, the aforementioned theoretical 
intensities were compared with the observed intensities using 
a least-squares minimisation method as follows: 



2 

X = 



E 



lobs (^0 ipred (^i) 
Ipred 



(4) 



where \ 2 is the sum over the combination of n lines used 
for either N II or O II, as a function of T c and 7V C . I b s is 
the observed intensity for the transition A; and I pre d is the 
predicted intensity. For a blended feature, the observed in- 
tensity is compared with the combined theoretical intensities 
of the blended lines, described in Eqn. (3). As an example, 
Figs. 6 and 7 show the log \ 2 distributions for 4 PNe (Hf 2-2, 
M 1-42, NGC 6153 and NGC 7009) for N n and O n, respec- 
tively. The location of the minimum log \ 2 value can be 
determined for each comparison, thus providing the optimal 
(log T c ,log 2V e ), as shown by the white cross hairs in both 
figures. 

2-4-2 Error estimate 

For many objects, the resulting temperatures and densities 
are well confined. However, it is difficult to decide how reli- 
able these diagnostic results could be. Not all of the surveys 
or individual detailed studies of PNe and/or H II regions pro- 
vided errors along with their observational intensities. This 
makes estimating the uncertainties for the optimal (log T c , 
log A c ) rather complicated. A few nebulae had observational 
errors presented in the literature and are adopted in the 
current analyses. For those nebulae whose measurement un- 
certainties are not presented in the literature, we estimated 
errors according to their observational conditions and data 
quality. In order to determine the uncertainties for the op- 
timal (log T c , log N e ), we carried out an extra set of calcu- 
lations that begin with random number generation. 

For each line in a combination of N II or O II ORLs (e.g. 
V3 A5666, V3 A5679, V39a A4035, and V39b A4041 for N n, 
and VI A4662, VI A4649, V48a A4089, and V48c A4087 for 
O n) that are used for plasma diagnostics, we used idl to 
generate a number, A a i m , of intensity values, J s im(A;), with 
a normal (Gaussian) distribution centred around the origi- 
nal observational intensity, 7 bs(Ai), along with a standard 
deviation obtained from the observational error, cr(/obs)> °f 
this distribution. 

The mean intensities of the simulations and standard 
deviations, I^ an and a(I™£ an ), are calculated as follows: 



iVeim 



and 



/ rmean\ 

<n4im ) = 



TtT (4 m (A) - (J sim )) 2 

Asim 



(5) 



(6) 



The standard deviation of these generated numbers also 
agree well with the observational errors. Each combination 
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a) loge(N 1 1 V3 X5666) + 26 




log e (X) + 26 b) log e (N 1 1 V3 X5679) + 26 
3.50 



log e {X) + 26 
3.581 




log e {X) + 26 
3.26 



2 3 
d) log 8 (N 1 1 V39a X4035) + 26 



4 

logN e 



log e {X) + 26 
3.121 




Figure 4. Emissivitics of the N II a) V3 A5666, b) V3 A5679, c) V39b A4041 and d) V39a A4035 lines as a function of electron 
temperature and density, as derived by Eqn. (1). The colour bar indicates the log e(A) + 26 values in units of ergscm+ 3 s" 1 . 



of the randomly generated line intensities, lsj m (A), where 
i is from 1 to N s im, is used to determine a unique opti- 
mal (log T , log N e ), following the method described in Sec- 
tion 2.4.1. 

Errors of the optimal log T e and log N e of each neb- 
ula are then calculated from these -/V s i m temperatures and 
densities, which can be denoted as Tg and Nl, where i = 1, 
7Vsi m . For those objects whose log T c - log N c diagram 
shows multiple peaks, standard deviations of T e and -/V e are 
calculated, which are centred on the mean values of these 
iV s i m temperatures and densities. The random intensities of 
each ORL are then normally distributed within the observa- 
tional errors of the observed intensities. The program would 
then calculate the log \ 2 along the entire T e and iV e grid 
for each combination of the random intensities and deter- 
mine the minimum value, thus deriving the unique optimal 
(log T c ,log N e ) for that simulation. 

Figs. 8 and 9 display the frequencies of these optimal 
log T and log iV c distributions for the N II and O II ORLs, 
respectively, for the same 4 PNe shown in Figs. 6 and 7. 
The errors for the optimal log T c and log N c are derived as 
follows: 



and 



1 v sin 



E -Wain, rpi 
i=l 1 » 

Nsim 



iVsim 



(7) 



(8) 



with respective standard deviations calculated as follows: 



/rr-imean \ 
"Vsim ) = 

and 

f-(iVsi m ) = 



iVsim 



(9) 



(10) 



These values are represented in Figs. 8 and 9 as green er- 
ror bars centred at T^ an and N Al ^ an , but projected on 
the white crosshairs of the optimal log T c and log iV loca- 
tions obtained from observations, thus providing upper and 
lower limits. The white solid contour designates the 1-a level 
around the peak. We project the standard deviations onto 
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. The same as Figure 4, but for O II a) VI A4649, b) VI A4662, c) V48a A4089 and d) V48c A4087 lines. 



the optimal log T c and log N c locations obtained from ob- 
servations because they are calculated using the l-cr errors 
of the observed intensities. 



3 RESULTS AND DISCUSSION 

3.1 Temperature structure of photoionized 
nebulae 

Figs. 10 and 11 show the electron temperatures from the lit- 
erature and the current ORL diagnostics for 46 PNe and 
for 42 H II regions, respectively. The PNe in Fig. 10 were 
chosen out of the whole list because they met the following 
criteria: 1) T e ([0 ill], T e (H i) and T G (He i) are available from 
literature, 2) 2.6 < log T c (0 II ORLs)< 4.2 and/or 3) 2.1 < 
log T e (N II ORLs)< 4.3. For the PNe plotted in Fig. 10, the 
mean value of log T e ([0 ill]) is 4.01±0.09 [K] (red solid line). 
The mean value of log T e (H I BJ) is 3.91±0.29 [K] (orange 
dotted line). The mean value of log T c (He I A7281/A6678) is 
3.67±0.19 [K] (purple dashed line). The mean value of 
log T c (0 II ORLs) is 3.22±0.52 [K] (green dot-dashed 
line). The mean value of log T C (N II ORLs) is 3.20±0.59 
[K] (blue dot-dot-dot-dashed line). In the case of the PNe, 



the mean values clearly show the temperature sequence: 
Te(ORLs) sC T c (He i) sC T C (H i) sC T c (CELs), which is 
consistent with predictio ns from the bi-abundance model 
jLiu et al.ll2000l : |Liul|2003h . While the physical conditions of 
the main nebula and cold, H-deficient component may vary 
from nebula to nebula, over all the nebulae there is a gen- 
eral trend that many contain a secondary cold, H-deficient 
component. 



For the H II regions plotted in Fig. 11, the mean value 
of log T c ([0 ill]) is 4.00±0.11 [K] (red solid line). The 
mean value of log T C (H I BJ) is 3.87±0.08 [K] (orange dot- 
ted line). The mean value of log T c (He I A7281/A6678) is 
3.91 [K] (purple dashed line). The mean value of 
log T c (0 ii ORLs) is 3.24±0.57 [K] (green dot-dashed 
line). The mean value of log T C (N II ORLs) is 3.15±0.87 
[K] (blue dot-dot-dot-dashed line). In the case of the H II re- 
gions, the mean values show the temperature relation: 
T c (ORLs) < T C (H i) < Te(Hei) < T c (CELs). The difference 
in this relation are most likely due to the faint nature and 
large measurement errors of He I lines in some of the H II re- 
gions. 
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7 0bs (^5666) + / obs (^5679) + 7 obs (A4041) + 7 obs (A4035) 



b) M 1-42 2 adf (O 2 l = 22.0: 

7 0bs (W5666) +X bs (k5679) + / obs (U041) + 7 obs (A4035) 




c) NGC6153 1 adf (O 2 *) = 9.2: 

7 0bs (^5666) + Q[X5679) + I obs (U041) + J obs (M035) 



logx 2 
1.10 



4 

logN, 

d) NGC 7009 21 adf (O 2 *) = 4.7: 

7 0bs (W5666) + / obs (W5679) + 7 obs (A4041) + 7 obs (A4035) 



logx 2 



1.01 




0.74 


m , 


0.47 




0.19 




-0.08 




-0.35 




-0.63 




-0.90 




-1.17 




-1.44 




-1.72 






23456 23456 
log N e log N e 

Figure 6. Logx 2 (T e ,N e ) distributions for N II V3 A5666, V3 A5679, V39b A4041 and V39a A4035 over the entire log T c - log N e grid 
for 4 PNe. The white crosshair pinpoints the absolute minimum x 2 value for the N II line comb inations in dicating the opti mal log T e 
and l og N e . The colou r bars indicate t he log \ 2 values. ADF values are f rom li terature: a) Hf2-2 1 llLiu et al.ll2006t) . b) M1-42 2 {l iuetajj 
120011) , c) NGC6153 15 llLiu et alfe oOO), and d) NGC 7009 21 l lFang fc LiulbOllT) . 



3.2 Physical evolution of photoionized nebulae 

Fig. 12 shows histograms of T c ([0 in]), T C (H I BJ), T c (He i), 
T c (0 n) and T C (N n), for PNe (blue) and H n regions (red). 
The He I ORL temperatures are derived from the line ra- 
tio A72 81/A6678 using the fitting functions of IZhang et al.l 
(2005b). Given the weakness of the N n and O n ORLs, the 
large scatter seen in the T c (0 n) and T C (N n) distributions 
are most likely due to observational uncertainties. The dia- 
gram clearly shows the relation of the average temperatures 
for the PNe, T e (ORLs) sC T c (He i) sC T C (H i) sC T c (CELs), 
which is also seen in Fig. 10. The average T e 's and stan- 
dard deviations are labelled for PNe (in red) and H n re- 
gions (in blue). Many PNe and H n regions in our analysis 
show T c (0 n) lower than 2.6 and T e (N n) lower than 2.1, 
the lowest T c 's for which the effective recombination coeffi- 
cients are available, respectively, although only a handful of 
H II regions had N II ORLs listed in their literature for us to 
perform our diagnostics. As can be seen in Fig. 11, the tem- 
perature histograms of [O III] CELs, He I, H I BJ, O II and 



N II ORLs show an increase in the width of the distribu- 
tions from [O III] through N II. This can be either caused by 
the observational uncertainties arising from increasingly dif- 
ficult measurements, especially for the extremely faint ORLs 
or possibly due to variations of physical properties of the 
cold H-deficient inclusions and main nebula from object to 
object, causing the scatter in T C (N n) and T e (0 n). 

Fig. 13 shows the ADF (0 2+ /H+) versus T e ([0 ill]) - 
T C (H I BJ), for PNe (blue) and H II regions (red). The 
123 PNe analysed here are shown as the blue open circles; 
amongst them the 46 PNe plotted in Fig. 10 are shown as 
the blue filled circles. A linear least-squares fit to the 46 
PNe from Fig. 10 with reliable electron temperatures indi- 
cates a positive correlation. The PNe shown in open circles 
and H II regions are not considered for the linear fit. This 
most likely indicates a tight correlation between the ADF of 
a nebulae and the temperature discrepancy. The 42 H II re- 
gions analysed here are shown as the red open squares with 
40 H II regions shown as red closed squares if they had pos- 
itive N II or O II log T c results. The H II regions plotted in 
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Figure 7. The same as Fig. 6 but for O II VI A4649, VI A4662, V48a A4089 and V48c A4087 lines over the entire 
for the same 4 PNe. 
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Fig. 13 just show large scatter, making it difficult to deter- 
mine if there is a correlation between the ADF value and 
temperature discrepancy. 



The nebulae with lower ADF values might be relatively 
young in their evolution processes when the main nebulae 
are very bright, and therefore the cold, H-deficient compo- 
nents might be mixed with the main nebula, and the ef- 
fects of the cold, H-deficient components are insignificant. 
However, for the nebulae with high ADF values, the nebula 
might have expanded such that the main nebula has a rel- 
atively low surface brightness compared to the condensed, 
cold H-deficient condensations, thus making it easier to dis- 
tinguish the two based on the temperatures derived from 
CELs and those derived from ORLs. Therefore for most 
PNe, there may be an evolution from younger, more compact 
nebulae where the cold, H-deficient component is embedded 
and almost indistinguishable from the main bright nebula 
to older, more diffuse nebulae where the fainter main neb- 
ula distinctly envelopes the cold, H-deficient condensations. 



3.3 Comments on several archetypal planetary 
nebulae 

3.3.1 Hf2-2 

There are several well-studied planetary nebulae that have 
been observed for many decades. Hf2-2 is a well-known 
south ern PN, that ha s an unusually high ADF (0 2+ ) value 
of 84 (|Liu et alfeOOot ). They observed this PN with the ESO 
1.52-m telescope with three separate slit widths of 2, 4, and 
8 arcsec, using the 2-arcsec width for maximum spectral res- 
olution. Hf 2-2 is known to have a close b inary system wit h 
a very short orbital period (0.398 days. IL"utz et al.lll998l ). 
It has also shown a very peculiar nature and may be in- 
cluded in the 'born-again' PNe scenario, relating the partic- 
ularly large ADF to the phenomenon of novae, as stated in 
I Wesson et ail d2003h. From their 2-arcsec slit-width observa- 
tions, iLiu et alJ ( 2006J) derived the following electron tem- 
peratures: T e ([0 in]) = 8710, T e (He I A6678/A5876) = 1570, 
and T C (H i) = 933 K. From the plasma diagnostics based 
on the N n ORL lines of V3 A5666, V3 A5679, V39b A4041, 
V39a A4035, we derived an electron temperature of 398li 35 
K, as shown in the first panels of Fig. 5 for the \ 2 minimisa- 
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Figure 8. 2-dimensional frequencies of the optimal log T c and log N c locations calculated from 10 000 simulations within a l-cr normal 
(Gaussian) distribution for N II V3 A5666, V3 A5679, V39b A4041 and V39a A4035 over a portion of the logT c - log N c grid for the 
same 4 PNe as Figs 5 and 6, sorted by descending ADF. The white crosshair pinpoints the optimal log T c and log N c location obtained 
from the observed intensities. The l ight-green error bars also sho w the standard deviations as deri ved from Eqn. 6 . The ADF values are 
adopted from liter ature: a) Hf2-2 X llLiu et al.ll2006h , b) M 1-42 2 jLiu et alfeOOlh , c) NGC 6153 15 jLiu et al.ll200oT) , and d) NGC 7009 21 
llFang fc Liull201l[) . 



tion and Fig. 8 for the standard deviations. From the plasma 
diagnostics based on the O n ORLs of VI A4649, VI A4661, 
V48a A4089, V48c A4087, we derived an electron tempera- 
ture of 3160lg2o 

K, as shown in the first panels of Fig. 7 for 
the % 2 minimisation and Fig. 9 for the standard deviations. 



3.3.2 Ml-42 

Another PN is the Ga lactic bulge nebul a M 1-42, which 
was first d iscove red by [Minkowski l| 19461 ) . More recently, 
iLiu et all (|200ll ) observed Ml-42 along with M2-36 us- 
ing the ESO 1.52-m telesc ope. For M 1-42 they derived 
a fairly high ADF of 22. ILiu et al.1 (|200ll ) derived the 
following electron temperatures: T e ([0 in]) = 9120 K, 
T e (He i A7821/A6678) = 3790 K, and T C (H i) = 3980 K. 
From our plasma diagnostics based on the N II ORLs, V3 
A5666, V3 A5679, V39b A4041 and V39a A4035, we derived 
an electron temperature of 3980^jggQ K, as shown in the 



second panel of Fig. 5 for the \ 2 minimisation and Fig. 8 
for the standard deviations. But as one can see in the two 
figures, the electron density of 10 6 cm -3 only provides an 
upper limit for the diagnostics. From our plasma diagnos- 
tics based on the O n ORLs, VI A4649, VI A4661, V48a 
A4089 and V48c A4087, we derived an electron temperature 
of 1410_i5o K, as shown in the second panel of Fig. 7 for the 
X 2 minimisation and Fig. 9 for the standard deviations. 



3.3.3 NGC 6153 

Over the last several decades, detailed studies have been 
put forth on the well-known planetary nebula NGC 6153, 
a possibly su p er-me t al-rich nebula a nd first noted by 
iPottasch et"ail l|l984h . ILiu et al.l (|2000h derived a moder- 
ately high ADF of 9.2. Deduced from their electron den- 
sity distribution a nd the optical appearance of NGC 6153, 
lYuan et al.l l|201ll ) suggested it is most likely a bipolar neb- 
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Figure 9. The same as Figure 8 but with O II VI A4649, VI A4662, V48a A4089 and V48c A4087 lines for the same 4 PNe. 



ula probably with a central cavity and a density-enhanced 
waist, as viewed at a large angle to its polar axis. However, 
they also argue that a convincing physical model account- 
ing for the full range of behaviour NGC 6153 exhibits is still 
to be found. From the strong emission lines, they derived 
the following electron temperatures: T e ([0 in]) = 7240 K, 
T e (He i A7821/A6678) = 3260 K, and T C (H i) = 6030 K. 
Based on the N n ORLs of V3 A5666, V3 A5679, V39b A4041 
and V39a A4035, our plasma diagnostics derived an electron 
temperature of 2000±^ K, as shown in the third panel 
of Fig. 6 for the x 2 minimisation and Fig. 8 for the stan- 
dard deviations. Based on the O n ORLs of VI A4649, VI 
A4661, V48a A4089 and V48c A4087, our plasma diagnostics 
yields an electron temperature of I78OI230 K, as shown in 
the third panel of Fig. 7 for the \ 2 minimisation and Fig. 9 
for the standard deviations. 



3.3.4 NGC 7009 

NGC 7009 has had hundreds of publications over the last 
50 years or more. Also known as the Saturn Nebula, this 
large, double-ringed, high-surface-brightness nebula is par- 
ticularly well known for having an unusually rich and strong 



O 11 optical permitted lines ever since the early high- 
resolu tion photographic spectro scopy observations of IWvsel 
l|l942h and lAller fc Kalerl (|l964T ) , having observed more than 
100 O 11 permitted t ransitions dLiu et al.l Il995). A much 
more recent work by iFang fc Liu ( 2011 ) has made use of 
very deep CCD spectrum of the Saturn Nebula covering 
from 3040 to 11 000 A. They derived a reasonably high ADF 
of about 5. They identified over 1000 emission lines with 81 
per cent attributed to permitted lines and more than 200 
O 11 permitted lines. Their observations were made using 
the ESO 1.52 m telescope and the William Herschel Tele- 
scope (WHT), and securing all spectra with a long slit. 
Table 3 in their paper lists a very complete set of elec- 
tron temperatures and densities derived from numerous line 
ratios. They derived the following electron temperatures: 
T e ([0 in]) = 9810 K, T (He 1 A7821/A6678) = 5100 K, and 
T C (H 1) = 6420 K. From our plasma diagnostics based on 
the N 11 lines of V3 A5666, V3 A5679, V39b A4041 and V39a 
A4035, we derived an electron temperature of 1260^235^ 

K. 

as shown in the fourth panel of Fig. 6 for the \ 2 minimisa- 
tion and Fig. 8 for the standard deviations. Due to the data 
quality, our x 2 minimisation technique only provides an up- 
per limit to the electron density for N II ORLs. But judging 
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Figure 12. Distributions of T a derived from a) [O III] CELs; b) 
H I Balmer Jump; c) He I line ratio A7281/A6678 from litera- 
ture and our own calculations; d) O II ORL diagnostics; and e) 
N II ORL diagnostics for PNe (blue) and H II regions (red). For the 
distribution in each panel, the average T c and standard deviations 
are labelled for PNe (blue) and H II regions (red). 

from the 2-D distributions of the results derived from the 
randomly generated simulated intensities, a secondary peak 
occurs at an electron density of 1580 cm -3 and a tempera- 
ture of 2000 K. From our plasma diagnostics based on the 
O n lines of VI A4649, VI A4661, V48a A4089 and V48c 
A4087, we derived an electron temperature of 1580lJJg K, 
as shown in the fourth panel of Fig. 7 for the x 2 minimisation 
and Fig. 9 for the standard deviations. 



3.3.5 Abell58 & Abell30 

Another rather peculiar planetary nebula is Abell58, 
a "born-again" planetary nebula, known to contain 
an H-deficien t knot surrounding V605 Aql, according to 
IWesson et al.l (|2008l ). Hidden by a thick dusty torus, 
V605 Aql cannot be seen directly. However, from the sur- 
face abundances d erived for the s tar fro m observations of 
its scattered light, IClavton et all (|200tfl_ found it to be a 
typical Wolf-Rayet central star of PN. IWesson et al.l (2008) 
did not measure any N n ORLs in Abell58. From the 
plasma diagnostics based on the O n ORLs of VI A4649, 
VI A4661, V10 A4075 and V10 A4069.62 (blended with V10 
A4069.89), we derived an upper limit to the electron tem- 
per ature (~ 15 800 K). 

IWesson et al.l (|2003T ) have provided detailed optical 
analysis of two of the H-deficient knots (Jl & J3) in another 
"born-again" PN Abell30. This nebula consists of a rather 
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Figure 13. Log ADF (0 2 +/H+) versus T c ([0 III]) - T C (H I BJ) re- 
lation for PNe (blue) and H II regions (red). The blue open circles 
are all 123 PNe listed in Table 1. The blue closed circles are all 
the 46 PNe plotted in Fig. 10. The solid black line shows a least- 
squares fit to the 46 PNe plotted in Fig. 10. The red open squares 
are all 42 H II regions listed in Table 2. The red closed squares are 
H II regions with either log T e (0 II) or log T C (N II) greater than 
zero. 

large, 120 arcsec across, spherical shell of low surface bright- 
ness a nd s everal bright clump s, first discovered by Ijacobvl 
l| 19791 ) and lHazard et alj (|l980h . of material about 10 arcsec 
away from the central star. For IWesson et all l|200ot ). they 
observed Abell 30 with the 4.2 m WHT at the Observatorio 
del Roque de los Muchachos, on La Palma, Spain, with a 
0.82 arcsec slit width. These two knots have extremely high 
AD F values of 598 and 7 66 for Jl and J3, respectively. 

IWesson et ail <|2003h derived the following electron 
temperatures for knot Jl: T e ([0 hi]) = 20 8 K and 
T e (He I A5876/A4471) = 350 K. IWesson et aD (|2003l ) did 
not measure any N II ORLs in the knot Jl. From the 
plasma diagnostics based on the O II ORLs of VI A4649, 
VI A4661, V10 A4075 and V10 A4069.62 (blended with V10 
A4069.89), we derived an upper limit to the el e ctron temper- 
ature (~15 800 K). For knot J3. 1 Wesson et all (120031 ) derived 
the following electron temperatures: T o ([0 III]) = 17 960 K, 
T e (He I A6678/A4471) = 9 240 K. They didn't measure any 
N II ORLs in this knot either. Plasma diagnostics based on 
the O II ORLs yielded an upper limit to the electron tem- 
perature (~15 800 K). 



3.4 Results and discussion 

Tables 1 and 2 list the electron temperatures and densities 
for the 167 nebulae, including 123 PNe and 42 H II regions, 
respectively. These objects are presented in descending or- 
der of ADF (0 2+ /H+). Column 1 is the nebula name. Col- 
umn 2 contains the ADF values from literature. Column 3 
is the log T o ([0 III]) [K] value from literature. Column 4 is 
the log T C (BJ) [K] value from literature. Column 5 is the 
log T e (N n) [K] result from our plasma diagnostics. Col- 
umn 6 is the log A r e (N n) [cm -3 ] result from our plasma 
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diagnostics. Column 7 is the log T o (0 n) [K] result from our 
plasma diagnostics. Column 8 is the log N c (0 n) [cm -3 ] re- 
sult from our plasma diagnostics. 

Detailed modelling for the He I recombination 
lines have been computed for a c ouple decades now. 
Benjamin. Skillman fc Smitsl jl999) combined the 
Smitsl fimxh models with collisional transitions of 
Sawev fc Berringtonl (11993}) producin g accura te emis- 
sivities. More recently^ [zhang et alJ l|2005bl ) applied 
computational models for the He I recombination lines 
for several dozen PNe. They also provide fitting functions 
spanning from 5 000 to 20 000 K with fitting parameters 
that are density-dependent. 



Tables 3 and 4 show the electron temperatures derived 
from the He I ORLs for PNe and H n regions, respec- 
tively. The three He I line ratios A7281/A6678, A6678/A4471 
and A6678/A5876 are used, and the objects are in descend- 
ing order of the ADF v alues. The PN diagnostic results 
from IZhang et al l l|2005bh are adopte d. For those P Ne tha t 
are not included in the sample of IZhang et al l (|2005bh . 
we derived elec t ron tem peratures using the fitting formula 
of IZhang et all (|2005bT l. and the He I line intensities are 
adopted from literature. Column 1 is the nebula name. Col- 
umn 2 is the ADF value from literature. Columns 3 through 
5 are the T c 's derived from the He I line ratios. Most litera- 
tures provided line intensities for the four He I lines A7281, 
A6678, A5876 and A4471. 

Resonance fluores cence may also have a minor effect on 
the strengths of ORLs. iGrandil fl976r) discussed in detail the 
effects of resonance fluorescence on permitted transitions. It 
has been known that the N II permitted lines from the low- 
lying 3d-3p and 3p-3s triplet arrays, whose upper levels 
are linked to the ground term 2p 2 3 P by resonance line s, can 
be enhanced by fluorescence excitation. IGrandil (jl976T j used 
photoionization models to study the excitation mechanisms 
of permitted transitions from common heavy element ions 
observed in the spectra of the Orion nebula and two Galactic 
PNe NGC 7027 and NGC 7662. 

IGrandil j 19761 ) found that while the N n M28 3d 3 D° - 
3p 3 P multiplet is excited by both recombination and 
continuum fluorescence of the starlight, emission of the 
N n M3 3p 3 D-3s 3 P°, M5 3p 3 P-3s 3 P° and M30 
4s 3 P°-3p 3 P multiplets are dominated by fluorescence ex- 
citation of the N II 4s 3 P° level by the He I ls8p 1 p°- 
ls 2 1 S A508.643 resonance line, which coincides in wave- 
length with the N II 2p4s 3 P? -2p 2 3 P A508.668 line. Flu- 
orescence excitation cannot excite the singlet transitions or 
transitions from the 3d-4f configuration. Given the fact that 
measurements of the N II ORLs for the majority of nebulae 
samples are of relatively large uncertainty due to weakness 
and blend, the effects of fluorescence on the plasma diag- 
nostic results are insignificant. For those nebulae with good 
observations, e.g., the several archetypal PNe discussed in 
Section 3.3, detailed modelling is needed to e stimate t he ex - 
act contribution of fluorescence mechanisms. IGrandil (|l976l ) 
showed that the dominant excitation mechanism of the 
strongest O II permitted lines is recombination. 



4 SUMMARY 

For nearly four decades, analyses of gaseous nebulae 
are mainly based on the bright CELs, whose emissiv- 
ities are acutely sensitive to electron temperature (e.g. 
lOsterbrock fc Ferlandl 120061 ). Thus the abundances derived 
from CELs could be subject to the uncertainties caused 
by temperature fluctuations if present in the nebulae. The 
ORLs of heavy element ions, although much weaker by na- 
ture, are much less sensitive to temperature, and the resul- 
tant abundances are supposed to be more reliable, provided 
that very accurate measurements of the heavy element ORLs 
can be obtained. In the current paper, we demonstrate a new 
plasma diagnostic method based on the N II and O II optical 
recombination spectra, using the new effective recombina- 
tion coefficients. This is the first work devoted to nebular 
analysis using heavy element ORLs, with a large sample of 
PNe and H II regions considered. The results show system- 
atic differences between the electron temperatures derived 
from CELs, the optical recombination spectra of H I and 
He I, and the N II and O II ORLs. The observed tempera- 
ture sequence has been found in previous observations, and 
is in agreement with the expe ctation of the bi-abundance 
nebular model (|Liu et al.ll2000i ). Although very deep, high- 
resolution spectra of gaseous nebulae, especially accurate 
measurements of heavy element ORLs, are still rare, we need 
to develop plasma diagnostic tools based on the heavy ele- 
ment ORLs for future study, given that the most compre- 
hensive treatment of the N II and O II recombination under 
the physical conditions of gaseous nebulae (i.e., the effective 
recombination coefficients for the N II and O II recombina- 
tion spectra) are now available. The main purpose of this 
paper is to demonstrate the new method of nebular analysis 
and show its potential application. 

The method presented in the current work for deter- 
mining T e and N e from multiple transitions of N II and 
O II is general and quite promising for the future deep spec- 
troscopic studies of gaseous nebulae. Transitions from mul- 
tiplets V3 and V39 for N II and from VI, V10 and V48 
for O II tend to be the strongest and most reliable to con- 
strain the x 2 distribution over the entire log T c - log N e grid. 
However, some other N II and O II transitions can be used 
as well. Therefore, this method can potentially cover broad 
wavelength ranges of a spectrum instead of just determin- 
ing the T c and N e from just a couple of lines, depending on 
the data quality. Given that the heavy element ORLs are 
intrinsically faint, high signal-to-noise ratios are needed to 
improve the data quality. 
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Table 1. Electron temperatures and densities for PNe. 
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^• zu 0.64 


3.80 ; 40 


NGC 6879 13 




2 


.46 


4.02 13 


3.93 13 






4.20 


z -'°0.02 


IC 4191 fixed 4 




2 


,40 


4.03 4 


4.02 4 


2.10 


5.90 


2.60 


2.00 


IC 4191 scanning 


A 


2 


,40 


4.03 4 


4.02 4 


2.10 


5.90 


2.60 


2 00 01 

z - uu o.oi 


NGC 3132 10 




2 


.40 


3.98 10 


4.03 10 


3 30 12 


6 OO - 49 


3 20 ' 34 
°- zu n .12 




NGC 6720 7 




2 


.40 


4.03 7 


3.90 7 


z - lu 0.03 


5 70 - 31 
°-' u 0.37 


2 fi=5°-° s 


2 20°' 07 
z - zu 0.06 



I iLiu et al.l ||2000|) ; 2 iLiu et al.l d200lh ; 3 iRuiz et all l l2003h ; 4 lTsamis et~al] d2003bh ; 5 IZhang fc Liul i2003Tl ; 

6 lErcolano etal] (120041 1; 7 ILiu et al.l | |2004[) ; 8 |Peimbert et al.1 J2004l1 ; 9 ISharpee et alj d2004h ; 10 iTsamis et alj ( l2004h ; 

II IWesson fc Liul ||2004T ); 12 iRobertson-Tessi fc Garnettl j2005T) ; 13 IWesson et all <2005l) ; 14 IZhane et af] 11200531) ; 

15 ILiu et al.1 ^OOrj l; 16 ISharpee et al.1 hoof ); 17 l\Vang fc Liul | |2007|) ; 18 iGarci'a-Roias et al.1 feOOSft ; 19 lQtsuka et al.1 i2010h ; 
20 IWesson et all i2003h ; 21 iFang fc Liul ^OllT) ; 22 lEsteban et al.1 J2002h ; 23 lTsamis et al.1 ^OOSal ); 24 lEsteban et al.1 \2004 ); 
25 IGarci'a-Roias et al.1 Hfloi); 26 IGarci'a-Roias et al.1 l|2005l l: 27 IGarci'a-Roias et al.1 l|2006h ; 28 IGarci'a-Roias et al.1 l|2007h ; 
29 lLopez-Sanchez et al.1 l|2007l ); 30 lEsteban etaD l|2009h ; 31 IWesson et all j200Sft : 32 ITsamis et aD (|2008h . 
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Table 1. (continued) 

From Literature From Current Work 

Object ADF (0 2 +)log T ([0 in]) log T c (BJ)log T C (N n)log T C (N n)log T C (Q n)log T c (0 II) 









[K] 


[K] 


[K] 


[cm 3 ] 


[K] 


[cm 3 ] 


NGC 6818 10 


2 


.40 


4.12 10 


4.08 10 


3.20 52 


6.00°; 37 


2 95 11 
z - 9o 0.12 


4 75O.09 
*• ,0 0.15 


IC 4191 nebula 10 


2 


.40 


4.03 10 


4.02 10 


3- SU 0.03 


2-70 ;i 2 


3.80 ; 07 


3 90°- 15 
J - yu 0.14 


IC 4191 fixed 10 


2 


.40 


4.03 10 


4.02 10 


q cnO.08 
a ' ou n .04 


2.70 - 1 ! 


q okO.05 
QQ3 


4 n5 - 36 

*' UJ .24 


IC 4191 16 


2 


.40 


4.00 16 


3.90 16 


2 70 ' 09 
05 


2 20° : ° 2 
16 


3.60g°f 


5.00 


NGC 3918 4 


2 


.30 


4.10 4 


4.09 4 


2 10° :o1 

fi-SJ 


5 90° : ° 4 
fi-fi 4 
5 90 21 


2.60°;o 4 


2.650; 12 


NGC 6884 7 


2 


.30 


4.04 7 


4.06 7 


u 0.03 


2.60 


2.40g : g| 


IC 5217 13 


2 


.26 


4.05 13 


4.08 13 






3 50 - 08 


s.oogil 


NGC 3242 10 


2 


.20 


4.07 10 


4.01 10 


4 00 11 

^• uu o.oi 


2 50° 11 
z - ou 0.06 


q qcO.05 
o,JJ 0.06 


5.oog : j» 


NGC 5882 4 


2 


.20 


3.97 4 


3.89 4 






2.60 


q in018 
O,lu 0.28 


M 1-74 13 


2 


.14 


4.01 13 


3.89 13 






q 7t;0.30 


4 4^0.53 


NGC 5882 32 


2 


.11 










2 95°- 72 
33 


2 40 1 94 
z, ^ u 47 


NGC 5882 10 


2 


.10 


3.97 10 


3.89 10 


3 90 07 


3 OO 60 


q ecrO.09 
J.OOq 2 




Me 2-2 13 


2 


.10 


4.04 13 


4.04 13 


Q SO 03 


2 20 ' 05 
z - zu 0.19 


4 OO 03 


5.oo°: 23 


NGC 7662 7 


2 


.00 


4.13 7 


4.09 7 


2.10 


5.90 


2.60 


2.00 


NGC 5315 10 


2 


.00 


3.95 10 


3.93 10 


3 - 80 0i44 


4.70 1 ; 03 


4 OO - 02 


4 qn0.35 
i.oUq 27 


NGC 6807 13 


2 


.00 


4.04 13 


4.00 13 






4 OO - 04 
uu 04 


q i cO.43 
° °0 31 


NGC 5307 3 


1 


.95 




4.03 3 






q i rfl.15 
° J 15 


2 70 60 
' u 34 


IC 4406 4 


1 


.90 


4.00 4 


3.97 4 


2.10 


5.90 


2.60 


2-2og:i? 

2-80° : ? 7 


IC 4406 10 


1 


.90 


4.00 10 


3.97 10 


2.10 01 


5 qn0 05 

06 


3.00°;°* 


NGC 6741 7 


1 


.90 


3.99 7 


4.15 7 


4 9n o:34 

*- zu 0.45 


2 on° :s4 


2.60 


2.20 


NGC 6826 7 


1 


.90 


3.97 7 


3.94 7 


2 30 25 
z ' ou 0.13 


5.70 ' 41 


2 90"' 01 
z - Ju 0.02 


2.00°;°1 


My Cn 18 4 


1 


.80 


3.86 4 








2.60 


2-15 .is 


My Cn 18 10 


1 


.80 


4.10 10 


4.09 10 


4.30 


2.00 


q ocO.03 
°- zo 0.03 


2 c;n 02 


NGC 3918 10 


1 


.80 


4 OO 0,21 


2 70 ' 28 
z -' u 0.39 


2.60 


2 OO - 09 
z ' uu 0.07 


NGC 7027 12 


1 


.80 


4.15 12 


4.08 12 






2.65 05 


2.00 


NGC 7027 14 


1 


.80 


4.10 14 


4.08 14 


3.90 ; 14 


3 10 143 
' J - lu 0.65 


4 OO - 03 

Q.1Q 


5.00 


NGC 5315 s 


1 


.74 


3.95 s 


3.93 s 


4 -30 . 10 


6.00 


4 20°- og 
*' zu 0.13 


4 qc;0.09 
^' a ' J 0.30 


NGC 6790 7 


1 


.70 


4.11 7 


4.15 7 


2.10 


5.90 


2.60 


2.20 


NGC 6790 12 


1 


.70 


4.11 12 


4.15 12 






2.60 


2.00 


Hu 1-2 7 


1 


.60 


4.29 7 


4.30 7 


2.10 


5.90 


2.60 


2.00 


NGC 6572 7 


1 


.60 


4.03 7 


4.01 7 


2 tin - 59 
z - ou 0.06 


5 40 0.29 


2 70 - 01 
Q-Q5 


3.35 007 


NGC 6572 12 


1 


.60 


4.01 12 


4.01 12 




2 60°- & 
Z - Du 0.01 


ac;0.43 
z ' oo 0.07 


NGC 6572 S 12 


1 


.60 


4.01 12 


4.01 12 






2.60 162 


2.30 


NGC 6891 13 


1 


.52 


3.97 13 


3.77 13 


2 10 - 01 
u 01 




4 20 - 10 


3.300 : 19 


NGC 5315 4 


1 


.40 


3.95 4 


3.93 4 


5.90 


2 90 - 10 

06 


5.00° : 93 


M 3-32 17 


1 


.15 


3.95 17 


3.65 17 


3 10 ' 20 


5 30 ' 47 


3 50 - 06 
J - ou 0.02 


3.35°:j:J 


M 3-33 17 


1 


.10 


4.02 17 


3.77 17 








q 9CO.O8 


IC 4699 17 


1 


.09 


4.07 17 


4.08 17 


2 60 1 11 
uu 23 


3 40 1 41 

93 


3-40 010 




NGC 6439 17 


1 


.09 


4.02 17 


4.00 17 


4 30 01 
01 


5.70° :i f 


4.00H 1 


3-55q ' -^2 


H 1-41 17 


1 


.08 


3.99 17 


3.65 17 


2 90° :95 
z - au 0.09 


4.80? : « 


q 4^0.04 
J -^ o 0.05 


2 10° :21 
z,lu 0.13 


M 3-7 17 


1 


.07 


3.88 17 


3.84 17 


4.30 


2.00 


q onO.24 


9 85 - 27 
z,oo 0.26 


NGC 6620 17 


1 


.06 


3.98 17 


4.00 17 


4 20 (M14 
*- zu 0.52 


3 00 05 


q qc;0.08 
°- oo 0.02 


q 7K0.19 


H 1-50 17 


1 


.05 


4.04 17 


4.10 17 






q i r0.02 


z ' oo 0.07 


H 1-54 17 


1 


.05 


3.98 17 


4.10 17 


4.30 


2.00 


q okO.32 
o,ZJ 0.20 


q 7K0.32 
J '' ,J 0.24 


M 1-29 17 


1 


.05 


4.03 17 


4.00 17 


o no 80 


3 40 1 ' 45 


2.60 


2.00 


M 3-21 17 


1 


.05 


3.99 17 


4.04 17 


4 30 02 
*- ou 0.02 


3 - 80 o:36 


2.60 


2.00 


H 1-35 17 


1 


.04 


3.96 17 


4.00 17 


4.30 


2.00 


q 2 n0.02 
o,zu 0.01 


4 75O.IO 
4 - '°0.05 


H 1-42 17 


1 


.04 


3.99 17 


4.00 17 


4.30 


2.00 


3.600;g 3 


q 7^0.43 
°- ,d 0.30 


M 2-6 17 


1 


.04 


4.00 17 


4.07 17 






3.800;? 3 


5.00 


M 2-27 17 


1 


.04 


4.08 17 


4.15 17 


4.30 


2.00 


q qn0.29 
o.ou 17 


2 90 - 18 
z,yu 0.20 


M 2-33 17 


1 


.04 


3.91 17 


3.85 17 


3 30 53 
°- ou 0.44 


2-5o :i 


2 60°- 02 
Z - Du 0.01 


J - DU 0.22 


M 2-42 17 


1 


.04 


3.93 17 


4.15 17 


4.20 ln 


4 -80l.88 


2.60 


2.00 


M 3-29 17 


1 


.04 


3.96 17 


4.03 17 


3.10 015 


5.30 32 


2.60 


2.00 


NGC 6567 17 


1 


.04 


4.02 17 


4.08 17 


4.30 


2.00 


2 80°- 1B 
z ' ou 0.04 


q nt;0.16 
O.UOq 22 



I iLiu et al.l ||2000|) ; 2 iLiu et al.l ll200lh ; 3 iRuiz et all | |2003|) ; 4 lTsamis et"all l l2003bh ; 5 IZhang fc Liul J2003h ; 

6 lErcolano etal] ||2004| ); 7 ILiu et al.l l2004) ; s |Peimbert et al.l d2004ll ; 9 ISharpee et alj \2004) : 10 iTsamis et alj \2004) : 

II IWesson fc Liul ||2004| ); 12 iRobertson-Tessi fc Garnettl d2005h ; 13 IWesson et all d2005h ; 14 IZhang et af] ll2005all ; 

15 ILiu et alj ||2006| ); 16 ISharpee et al.l hoof ); 17 l\Vang fc Liul | |2007|) ; 18 iGarci'a-Roias et alj feOOSl) ; 19 lQtsuka et alj J2010l ); 
20 IWesson et al] (120031 ); 21 iFang fc Liul l l201ll) ; 22 lEsteban et al.l | |2002|) ; 23 lTsamis et al.l l l2003al ); 24 lEsteban et alj (120041 ); 
25 IGarci'a-Roias et alj | |2004| ); 26 IGarci'a-Roias et alj ||2005| ); 27 lGarci'a-Roias et alj J2006T) ; 28 IGarci'a-Roias et alj (l2007f) ; 
29 lLopez-Sanchez et alj d2007t) ; 30 lEsteban etal] d2009h ; 31 IWesson et al] ||2008|) ; 32 ITsamis et al] ll2008h . 
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Table 1. (continued) 



From Literature 



From Current Work 



Object 



ADF (0 2 +) log T e ([0 in]) log Te(BJ) log T C (N n) log T C (N n) log T c (0 n) log T e (0 n) 







[KJ 


l K \ 


IC 4593 12 


1.03 


3.92 12 


3.92 12 


He 2-118 17 


1.03 


4.10 17 


4.26 17 


M 1-61 17 


1.03 


3.95 17 


3.98 17 


M 2-4 17 


1.03 


3.93 17 


3.90 17 


NGC 6565 17 


1.03 


4.01 17 


3.93 17 


Vy 2-1 17 


1.03 


3.90 17 


3.94 17 


IC 418 9 


1.02 


3.95 9 




Cn 1-5 17 


1.02 


3.94 17 


4.00 17 


Cn 2-1 17 


1.02 


4.01 17 


4.03 17 


M 1-20 17 


1.02 


3.99 17 


4.08 17 


M 2-23 17 


1.02 


4.08 17 


3.70 17 


IC 2501 16 


1.01 


3.98 16 


3.85 16 


M 2-39 17 


0.95 


3.91 17 


3.74 17 


Cn 3-1 13 




3.88 13 


3.71 13 


M 2-31 17 




3.99 17 


4.15 17 



[K] 



4.30 

4 no 0, 15 
4 no - 01 

4.30 
4.30 
4 20 01 

o qnO.09 
°- ou 0.07 

U 0.15 



3. 10 



4.30 
4 20° 03 
4 in - 06 

4.30 



2.00 

2 -30 ( , 76 
2.00 
2.00 
2.00 

6 00 42 
d.uu 54 

6.00 



2.00 

5.80 
2.00 



2 80 ' 17 
3.05 05 

3.00 

4.20 

q qcO.63 
4 20 S - 02 

4.20 

2.60 

q oc:0.15 

2 60°- 01 
4 150-06 

4 00°- 06 
2 R50.04 



2 00 - 47 

2 45 - 21 
9 okO.02 

Q-Q6 

3 450-04 

q fi50.33 

q 9 c;0.05 

q cc;0.28 
°- oo 0.16 

2.00 
q 7=0.78 
0.29 

4 OO 013 

5 00°- 05 
°' UU .94 

5 00°- 02 
°- UU 0.02 

2 60 05 



9 ISharpee et al.l i2004h : ^Wesson etafl j2005h : 1( lSharpee et~all <2007ft : 17 IWang fc Liul ||2007|) . 
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Table 2. Electron temperatures and densities for H II regions. 



From Literature From Current Work 

Object ADF (0 2 +) log T c ([0 m]) log T C (BJ) log T C (N n) log T C (N n) log T (0 n) log T (0 n) 

[K] [K] [K] [cm- 3 ] [K] [cm" 3 ] 



SMC NS7 1 


2 


.80 


4.09 1 








2.60 


2.20 


M 17 4 


2 


.66 


3.91 4 


3.89 4 


4-00 ,- 1 ;? 

0.41 


2.40 n ,„ 

0.79 


2.60 


2.00 


LMC N141 1 


2 


.60 


4.07 1 








2.60 


2.00 


NGC 3576 4 


2 


.21 


3.95 4 


3.91 4 


3.30 01 


0.03 


2.60 


2.00 


30 Doradus 4 


1 


.76 


4.00 4 




3.30°,-^ 


6.00S?i 

0.29 


2.60 


2.00 


LMC N66 1 


1 


.61 


4.26 1 












LMC Nil B 4 


1 


.47 


3.97 4 








2 60 01 
i,uu 0.01 


2 OO 02 
z - uu o.oi 


PB 6-2 2 


1 


.12 


4.18 2 












PB 6-1 2 


1 


.10 


4.20 2 












PB 8 2 


1 


.05 


3.84 2 


3.71 2 


0.13 


5.80 ,- 44 

0.32 


4.200? 2 

0.13 


3 30 12 
0.09 


NGC 2366 NGC 2363 3 


1 


.04 


4.20 3 








0. 14 


2.90S- S7 

0.55 


NGC 3603 s 


1 


.04 


3.96 s 








2.600-01 

0.01 


2 ^5 - 02 

0.09 


NGC 2867-2 2 


1 


.03 


4.06 2 


3.95 2 






2.60 


3 80°- 21 

O,ou 0.39 


S 311 7 


1 


.03 


3.95 7 


3.98 7 






2.60 


2.00 


M101 NGC 5461 3 


1 


.03 


3.93 3 








o ic0.35 
0.32 


4.65S-2 2 

0.76 


M101 NGC 5471 3 


1 


.03 


4.15 3 








2.85 - 95 


4.40, 7fi 


M 42 5 


1 


.02 


3.92 s 


3.90 5 


0.02 


2.300-? 9 

0.05 


4.20S!-°,g 
0.08 


3.550-1; 

0.14 


NGC 2867- 1 2 


1 


.02 


4.07 2 


3.95 2 






2.60 


3 . 75i i n 
1 .40 


NGC 5253 HII-2 10 


1 


.01 


4.08 10 








3 50S i? 


°- iU 0.33 


NGC 5253 UV-1 10 


1 


.01 


4.04 10 








4.200-0/; 

U. (0 


2.85,, cr 

O.oo 


NGC 5253 HII-1 10 


1 


.00 


4.08 10 








3 909.ll 

°- ayJ 0.67 


3 50 ,- 42 . 


NGC 5253 UV-2 10 


1 


.00 


4.04 10 












M33 NGC 604 3 


1 


.00 


3.91 3 








3 oo ,- , 3 . 


2 35O IO 


M33 NGC 604 11 


1 


.00 


3.91 11 








2.600-60 

O.zO 


2-ooS-ii 

U.09 


M101 NGC 5461 11 





.99 


3.93 11 








o qrr0.38 
0.24 


2.650- 2 } 

0.11 


NGC 2366 NGC 2363 11 





.97 


4.20 11 








^'^^0.68 


2.00S-?§ 

wv 0.12 


M101 NGC 5447 11 





.97 


3.93 11 


3.82 11 






3 05°- 22 
J.UOq 22 


9 1 c0.52 
z ' lo 0.22 


SMC N66 4 





.94 


4.09 4 








2.60 


2.00 


NGC 3576 6 





.91 


3.95 6 


3.91 6 


2.10 


5.90 


3 50°- 07 


9 7CO.20 

z - ,o 0.12 


NGC 4395 Reg 70 11 





.77 


4.03 11 








o qnO.16 


45O.17 
z -^°0.22 


NGC 2403 VS 38 11 





.61 


3.94 11 








4 05 002 
*- uo .32 


2 7 n0.02 


M 16 s 





,15 


3.88 s 


3.74 s 






3 OO 013 
J.UUg 27 


2 fid 1 - 00 
Z.DUq 41 


NGC 2403 VS 24 11 





,15 


3.91 11 








3 7n«- 3 
°- ,u 0.09 


2 4n°- 36 

z -^ u 0.19 


NGC 1741 Zone C 11 





.38 


3.93 11 








3.65 - 40 


2.05 - 69 


M 8 9 





.37 


3.91 9 


3.85 9 


2.10 


5.90 


■? i k0.33 
o.io 14 


Z.DUq 2 6 


M101 H1013 11 





.36 


3.87 11 








q cc0.44 


2 cqO.20 
Z.OUq 23 


M33 NGC 595 11 





.31 


3.87 11 








3 70° - 3a 
a - ,u 0.18 


2 ccO.34 
Z.OO 2g 


M 20 8 





.33 


3.89 s 


3.78 s 






2.60 


2.00 


NGC 2403 VS 44 11 





.30 


3.92 11 








q cc0.40 
O.OOq 29 


2 no - 18 

z - uu 0.09 


M 17 9 





.27 


3.91 9 




4.30 


2.00 


3 60°- 42 


2.300;| 2 


NGC 4861 BrightHII 11 





.27 


4.11 11 








4 -20 . 59 


2 -90 .64 


M31 K932 11 





.21 


3.97 11 








o 7 c0.36 
°- '°0.66 


2-000; 3 | 



1 iTsamis et all l|2003bl l: 2 IGarcfa-Roias et alj (200^ : 3 lEsteban et alj d2002t ); 4 iTsamis et al.l [|2003ah : 5 lEsteban et alj l|2004h ; 
6 iGarci'a-Roias et al.l ll2004h : 7 lGarci'a-Roias et al.l (|2005h : 8 IGarcfa-Roias et al.l l|2006l l: 9 iGarci'a-Roias et alj j2007h : 
1° lLopez-Sanchez et al.l d2007t) ; 11 lEsteban etafl J2009l> . 
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Table 3. The He I diagnostics from literature and current analysis for ^g^jg 3 (Continued) 
p N6i a e on lnue 



T e (He I) [K] 



T e (He 1) [K] 



Object 


ADF(0 2+ ) 


A7281 


A6678 


A6678 






A6678 


A4471 


A5876 


AK^ll QO 7^20 
PVDC11 oU Jo 


7ftft o 




yusu 


0Z0U 


A Qfi T1 20 

A Dell oU J 1 


£oq n 
Oyo.U 




000U 




I\ Dell Oo 


so. o 




^87 


1 9^0 

1Z0U 


hi z-z ^z circbec j 


°.a n 




1 070 
1U / U 


1 ^70 
10 / u 


TJf 9 9 //I ^ ,.^0^1 15 

hi z-z (^4 drcsccj 


c/i o 

o4.U 






QA O 
o4U 


TJf 99/Q „ v^co^l 15 

hi z-z ^0 drcsccj 


c/i o 

o4.U 




Q 1 ft 

olu 


1 /IftO 
14DU 


IN^jO 10U1 


Q9 o 
oz.U 




m onG 

01UU 


A T OOG 
41UU 


A/T 1 492 

1V1 ±-^iZ 


99 o 
zz.u 


^700 

/ yu 


9^00 

zouu 


ZOOU 


IN 4U 


17 Q 

1 / .0 




7/1 QO 
/ 4oU 


KA 70 
04 / U 


Ayr 9 9zi5 

IVl Z- Z4± 


1 / .U 






A^OO 


IN^Lj zUzz 


1 ft n 

10. u 




ly ( UU 




Mpp 909910 
IN ut' ZUZZ 


1 ft o 
10. u 




loyuu 




JJQJJlH 1 


11.0 






Qe;nol3 
OOUU 


IN OlOO 


9 

y.z 


OZDU 


1 0R00 
1UDUU 


^^70 
00 / U 


ir* 9ooq13 
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N II & O II ORL Diagnostics 



Table 4. Electron temperatures of H II regions derived from the He I recombination lines. 
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Object ADF(0 2+ ) A7281 A6678 A6678 

A6678 A4471 A5876 
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1.0 
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1.0 
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M 42 4 
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0.9 
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0.4 
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0.4 
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0.4 
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0.3 


2110 


6560 
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0.3 
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0.3 
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5590 


6690 
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